Measurement of long-range wavefunction correlations in an open microwave billiard 
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We investigate the statistical properties of wavefunctions in an open chaotic cavity. When the 
number of channels in the openings of the billiard is increased by varying the frequency, wavefunc- 
tions cross over from real to complex. The distribution of the phase rigidity, which characterizes 
the degree to which a wavefunction is complex, and long-range correlations of intensity and current 
density are studied as a function of the number of channels in the openings. All measured quantities 
are in perfect agreement with theoretical predictions. 
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From a statistical point of view, eigenvalues and eigen- 
functions of the wave equation in a chaotic billiard are 
well described by random matrix theory 1]: Depending 
on the presence or absence of time-reversal symmetry, 
they show the same distribution as eigenvalues and eigen- 
vectors of a large hermitian matrix with random Gaus- 
sian distributed real or complex elements. Whereas ran- 
dom matrix theory predicts a very characteristic eigen- 
value distribution with eigenvalue repulsion and spec- 
tral rigidity Q — , its predictions for eigenvectors are 
rather "uninteresting": eigenvector elements are inde- 
pendent Gaussian distributed real or complex random 
numbers. 

The situation is different for the eigenvectors of a ran- 
dom matrix that interpolates between the standard en- 
sembles with real and complex matrix elements. An ex- 
ample of such an interpolating ensemble is the "Pandey- 
Mehta" Hamiltonian 



H(a) =H + aHx, 



(1) 



where Hq and Hi are real and complex random hermi- 
tian matrices, respectively, and a is a crossover parame- 
ter. For such an interpolating random matrix ensemble 
the eigenvector elements have a non-Gaussian distribu- 
tion [j, and acquire correlations, both between ele- 
ments of the same eigenvector 0,0 and between different 
eigenvectors @- 

Experimental verification of this eigenvector distribu- 
tion has proven problematic because the long-range cor- 
relations and the deviations from Gaussian distributions 
are only of the order of a few percent 0, H, IE 0- 
An additional complication arises from the fact that the 
crossover parameter a needs to be fitted to the exper- 
iment. These complications may explain why measure- 
ments of wavefunction distributions in two-dimensional 
microwave cavities in which time-reversal symmetry was 
broken using magneto-optical effects were inconclusive 
with respect to the functional form of the probability 
distribution and did not reveal long-range wavefunction 
correlations 91. 



An alternative method to observe the real-to-complex 
crossover is to consider trav eling waves in a billiard that is 
opened to the outside world [lfj. For microwaves, such an 
open billiard is obtained by connecting a two-dimensional 
microwave cavity to waveguides. The parameter govern- 
ing the crossover from real to complex wavefunctions is 
the total number of channels N in the two waveguides. 
As was shown by one of the authors , wavefunctions 
in this crossover have a non-Gaussian distribution and 
long-range correlations, just like the eigenvectors of the 
Pandey-Mehta Hamiltonian JJJ. The main difference, 
however, is that the crossover parameter N is discrete 
and can be measured independently. This allows a fit- 
parameter free comparison of theory and experiment. 

We here report on the first measurement of such long- 
range wavefunction correlations in the real-to-complex 
crossover. The basic principles of the experiment can be 
found in Ref. We used a rounded rectangular cav- 
ity (2 f cm x 18 cm) coupled to two waveguides of width 
3 cm with a cut-off frequency at vt — 5 GHz. To break 
the symmetry of the shape of the resonator and to block 
direct transport, two half disks with a radius of 3 cm 
were placed in the resonator. To avoid unwanted reflec- 
tion, absorbers were placed at the end of the leads. We 
scanned the billiard on a square grid of 5 mm with a mov- 
able antenna A\ and measured transmission S\i in the 
range of 4- 18 GHz from a fixed antenna A2 in the end 
of the right lead. The fixed antenna had a metallic core 
of diameter 1 mm and a teflon coating while the probe 
antenna A\ was a thin wire of diameter 0.2 mm to mini- 
mize the leakage current. The lengths of the antenna A\ 
and antenna A2 were 4 and 5 mm respectively. 

For microwave frequencies v < v max — c/2d — 
18.75 GHz, where c is the velocity of light and d is the 
resonator height, the billiard is quasi- two-dimensional. 
In this regime there is an exact correspondence be- 
tween electrodynamics and quantum mechanics, where 
the component of the electric field perpendicular to the 
plane of the microwave billiard E z corresponds to the 
quantum- mechanical wave function tp. We normalize the 
wavefunction ip such that f dr\tfj(r)\ 2 — 1. Then the 
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FIG. 1: Left panel: Grey scale plot of measured intensity \ip\ 2 
at v — 12.015 GHz. Black corresponds to maximum inten- 
sity. Right panel: Measured current density for the same 
wave function . Arrow lengths indicate the magnitude of the 
Poynting vector. 



square of the electric field and the Poynting vector map 
to the normalized "intensity" and "current density" , re- 
spectively, 

7(r)=A|^(r)| 2 , j(r) = ^In#*(r)W(r)], (2) 

where A is the area of the billiard and k the wavenumber. 
Figure^shows typical intensity and current patterns thus 
obtained. 

The correspondence between the perpendicular com- 
ponent of the electric field E z and the wavefunction ip 
has been used previously to study the spatial distribu- 
tions and correlation functions of currents and vortices 
in open billiards [l3l ITU . In the present work we study 
spatial correlation functions of the squares of intensities 
and currents. The statistical average is taken by aver- 
aging over both position and frequency. It is for these 
quantities and for this full ensemble average that long- 
range correlations were predicted 0] . Before we describe 
the experimental results, we briefly summarize the con- 
clusions of Ref . Ull 

InRef.O the wavefunction distribution was described 
as the convolution of a Gaussian distribution with corre- 
lated real and imaginary parts and that of a single com- 
plex number p, the dot product of the wavefunction ip 
and its time reversed, 



p = dr -0(r) 



(3) 



The absolute value \p\ 2 is known as the "phase rigidity" 
of if; Q. The Gaussian wavefunction distribution at a 
fixed value of p implies a generalized Porter- Thomas dis- 
tribution for the intensity, 



P P (I) = 



: exp 







W 






U-\p\ 2 \ 



(4) 



so that the full intensity distribution is obtained by con- 
volution of Eq. (0J and the distribution p{p) of p, 



P(I)= J dpp(p)P p (I). 



(5) 



The distribution p(p) was calculated in Ref. I 111 using ran- 
dom matrix theory. 

In order to explain the origin of long-range correlations 
of intensity and current density in an open chaotic bil- 
liard, we consider the joint distribution of intensities at 
points r and r' with separation k\r — r'\ 3> 1, 



P[J(r),/(r')] = J dpp(p)P p [I(r)]P p [I(r')}. 



(6) 



For an open billiard, p has a nontrivial distribution, hence 
the long-range correlations of P[I(r), I(r')). Whereas 
random matrix theory predicts the long-range correla- 
tions of intensities through Eq. ©, it cannot alone pre- 
dict the long-range correlations of current densities and 
the precise dependence of these correlators on the separa- 
tion |r — r'| for fc|r — r'| of order unity. However, as shown 
in Ref. 111! the latter can be obtained from the random 
matrix result by making use of Berry's ansatz 15], which 
expresses tp as a random sum over plane waves, 



zk-r 



(7) 



Here the plane wave amplitudes a(k) have a Gaussian 
distribution with zero mean and with variance 



(a(k)o(-k)>=p(a(k)o*(k)), 



(8) 



where p is the (random) phase rigidity of ip. Performing 
the ensemble average using Eqs. (JJJ and |(HJl, correlators 
of intensity and current density are then expressed in 
terms of moments of the phase rigidity \p\ 2 . Long range 
correlations are found for correlators involving the square 
of the intensity and the current density, 



{i{vyi{v'Y) c 

{I{v) 2 J{v'f)c 
{J{vfj{T'f)c 



var|p| 2 + 4/ 2 (4 +13|p| 2 + |p| 4 } 
+ 4/ 4 (l + 4H 2 + |p| 4 ), 

-ivar|p| 2 + / 4 (2-H 2 -H 4 ), 
ivarH 2 + l/ 2 (l-2H 2 + H 4 ) 
+ J/ 4 (3-5|H 2 + 2|p| 4 ). (9) 



Here / = Jo(fc|r — r'|), J being the Bessel function, 
and the subscript "c" refers to the connected correlator, 
(AB) C = (AB) - (A){B). The relevant moments of the 
phase rigidity \p\ 2 were calculated in Ref. Hit For N = 2, 
4, and 6, one has {\p\ 2 ) = 0.7268, 0.5014, and 0.3918, and 
(|p| 4 ) = 0.6064, 0.3285, and 0.2155, respectively. 

We now describe the measured wavefunction distribu- 
tions. We first discuss wavefunction distributions mea- 
sured at a fixed frequency, and compare to the theory for 
the corresponding fixed value of p. A statistical distri- 
bution at a fixed frequency is obtained by varying the 
position of the antenna only. The phase rigidity \p\ 2 
can be measured independently using Eq. J2Jl. Fig. [21 
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FIG. 2: Intensity distribution for the wave function at v — 
12.015 GHz. Theoretical values from Eq. J1J are plotted with 
solid line. 

shows the intensity distribution for the intensity pattern 
shown in the left panel of Fig.ntogether with the theory 
of Eq. Q, using the measured value of \p\ 2 = 0.5202. 
As was discussed in Ref. there are frequency regimes 
where the leakage to the probe antenna becomes intolera- 
bly high, either due to the fact that the transport through 
the cavity is small, or due to some strongly scarred wave 
functions. In all such cases there were strong deviations 
from the generalized Porter-Thomas behavior described 
by Eq. (0J. We therefore only used frequency regimes 
where P(|V>| 2 ) was in agreement with theory on a confi- 
dence level of 90 percent. To have a well defined number 
of transversal modes within each waveguide, we investi- 
gate the frequency regimes 5-9.5, 10-14.5, 15 -18 GHz, 
corresponding to a total number of channels 7V=2,4, 6, 
respectively (i.e., 1, 2, and 3 propagating modes in each 
waveguide) 

We now describe our results for a full ensemble aver- 
age, in which both the position of the detector antenna 
and the frequency are varied. Fig. [3] shows the measured 
phase rigidity distribution P(|p| 2 ) together with the the- 
ory of Ref. |lli for N — 2, 4, and 6. Good agreement is 
found between experiment and theory, especially as there 
is no free parameter. Figures ^ and EH show mea- 
surement and theoretical prediction for the correlation 
functions of the squared intensity and the squared cur- 
rent density at positions r and r'. Since these correlation 
functions depend on the positions r and r' through the 
combination k\r— r'\ only, results from different frequency 
regimes can be superimposed by a proper scaling. In our 
analysis, we selected frequency windows for averaging of 
Av = 0.3 GHz guaranteeing that Av <C c/2vL, where L 
is the billiard size and c the velocity of wave propagation. 
The gaps in the N = 4, 6 histograms for small values of 
k\r — r'| reflect the spatial resolution limited to 5 mm due 
to the chosen grid size. For the long-range correlations, 
we observe excellent agreement between experiment and 
theory. The short-range oscillations predicted by theory 
are suppressed to a large extent in the experiment for 
the correlators that involve the current density, however. 
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FIG. 3: Phase rigidity distribution P(|p| 2 ) for different values 
of the total number of channels N. Solid curves indicate the 
theory of Ref. [U 
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FIG. 4: Connected correlator of squared intensity 
(/(r) 2 J(r') 2 ) c versus k\r — r'| for different values of the to- 
tal number of channels TV. 
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FIG. 5: Connected correlator of squared intensity and 
squared current density (7(r) 2 J(r') 2 ) c versus fc|r — r'| , for 
different values of the total number of channels N. 



This, again, is a consequence of the limited resolution. 
We would like to stress, however, that the asymptotic 
values of all correlations in the limit k\r — r'| — > oo is 
different from zero and in perfect agreement with the 
predictions from Eq. ©. 

Long-range correlations have also been observed in 
the transmission of microwaves 0] and visible light |J3 
through two- and three-dimensional disordered media, in 
which wave propagation is diffusive. There are important 
differences between the long-range correlations observed 
in Refs. 0, and those in ballistic systems, which are 
reported here. First, in Refs. 0, 0, long-range corre- 
lations appear already for the intensity autocorrelation 
function, whereas we find long-range correlations for cor- 
relators of second and higher moments only. Second, 
in Refs. flfil IrR the correlations scale with the inverse 
of the sample conductance, whereas in the present work 
there is no small parameter that sets the size of the long- 
range correlations, irrespective of sample size or distance, 
fn that sense, only the latter correlations are truly long 
range. 

This work was supported by the DFG, the NSF under 
Grant No. DMR 0334499 and by the Packard founda- 
tion (PWB). 




FIG. 6: Connected correlator of squared current density 
( J(r) 2 J(r') 2 ) c versus fc|r — r'| for different values of the total 
number of channels N. 
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